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MTrROVASClJI ^ BLOOD VOLUME MAGNETIC RESONANCE IMA CmG 

Background of the Invention 

The foUowing relates to the diagnostic imaging arts. It finds particular 
application in non-mvasive measurement by magnetic resonance imaging of cerebral blood 
volumes, and will be described with particular reference thereto. However, it also finds 
application in measurement by magnetic resonance imaging of blood volumes in other tissues. 

In-situ measurement of blood volume is usefiil in various clinical, diagnostic, 
and research appUcations. Local cerebral blood volume changes, for example, correlate with 
local neuronal activity in the brain. Cerebral blood volume measurements during physiological 
stimulation thus provides a tool for fonctional studies of brain activity. Cerebral blood volume 
measurements can also provide information about impaired and/or damaged tissue in stroke 
victims, as well as about lesions in many disorders, including, but not limited to cancer, 
vascular disorders, and the like. Blood volume imaging of other organs besides the brain can 
similarly provide functional and diagnostic data that is useful in clinical studies, diagnoses, and 
tests (e.g. stress tests or tests of vascular compliance). 

Various imaging modalities have been used to measure blood volume 
characteristics, including positron emission tomography, single photon emission computed 
tomography, and magnetic resonance imaging. Most of these are invasive approaches, in which 
the subject is administered a suitable contrast agent that selectively enhances blood contrast in 
the selected imaging modality. For magnetic resonance imaging, various paramagnetic contrast 
agents are commonly used for this purpose. The requirement of an administered contrast agent 
is a substantial disadvantage of these techniques. 

Magnetic resonance imaging of blood oxygenation level dependence (BOLD) 
is a non-invasive technique for indirectly measuring blood volume. In this technique, blood 
hemoglobin is used as an endogenous contrast agent. In one BOLD imaging approach, 
magnetic resonance imaging is performed as a fimction of physiological stimulation that causes 
changes in blood oxygenation level. Blood volume is estimated from BOLD measurements by 
making assumptions pertaining to other parameters that affect blood oxygenation level, such as 
blood flow. Hence, BOLD does not provide a direct measure of the blood volume. 

A disadvantage of both the invasive techniques and the BOLD techniques as 
applied to blood volume measurement is that these existing techniques generally do not 
differentiate between blood in large blood vessels, on the one hand, and perfused blood in 



small capiUaries or other microvessels, on the other hand. The blood volume in larger blood 
vessels is principally controlled by sympathetic regulation. In contrast, blood volume in 
microvessels having typical diameters of less than about 200 microns tends to vary to maintain 
local homeostasis or in response to chemicals such as vasodilators or vasorestrictive 
compounds. Consequently, the blood volume of microvessels responds to physiological 
perturbations such as local neuronal activity. For functional magnetic resonance imaging, the 
volume of blood in the microvessels is typically of principle interest, while the blood signal 
from larger blood vessels is interfering and thus undesirable. On the other hand, total 
microvascular plus macrovascular blood volume may change in some diseases, including but 
not limited to for instance arteriovenous malformations. 

The present mvention contranplates an improved apparatus and method that 
overcomes the aforementioned limitations and others. 

Brief Summary of the Invention 

According to one aspect, a magnetic resonance imaging method is provided. A 
blood-nulling magnetic resonance excitation sequence is performed that substantially nulls a 
magnetic resonance signal from blood. Subsequent to the performing of the blood-nulling 
magnetic resonance excitation sequence, a readout magnetic resonance sequence is performed 
to acquire a magnetic resonance signal from tissue other than the nulled blood. 

According to another aspect, a magnetic resonance system is disclosed. A 
blood nulling means is provided for performing a blood nulling magnetic resonance excitation 
sequence that substantially nulls a magnetic resonance signal from blood. A readout means is 
provided for performing a readout magnetic resonance sequence to acquire a magnetic 
resonance signal from tissue other than the nulled blood, the readout means operating 
subsequent to operation of the blood nulling means. 

One advantage resides in measuring the parenchymal vascular space 
occupancy, which is more sensitive to physiological perturbations tiian is the total vascular 
volume which includes the large blood vessels. 

Another advantage resides in providing measurements of the absolute blood 

voliune. 

Yet another advantage resides in providing images with substantially nulled 
blood magnetic resonance signal, i.e. a blood signal reduction sufficient to have MRI signal 
remaining that is predominantly from other tissues. 

Numerous additional advantages and benefits will become apparent to those of 
ordinary skill in tiie art upon reading the following detailed description of tiie preferred 
embodiments. 



Brief Description of the Drawings 

The invention may take form in various components and arrangements of 
components, and in various process operations and arrangements of process operations. The 
drawings are only for the purpose of iUustrating preferred embodiments and are not to be 
construed as limiting the invention. 

FIGURE 1 diagrammatically shows a magnetic resonance imaging system for 

measuring vascular space occupancy. 

FIGURE 2 diagrammatically shows a magnetic resonance imaging pulse 
sequence employing inversion recovery to null the blood signal. In FIGURE 2, the time 
intervals TI, TE, TR are not drawn to scale. 

FIGURE 3 shows a block diagram of a processor for computing a blood 
volume at rest and a blood volume change occurring responsive to a physiological perturbation. 

FIGURE 4 shows a block diagram of a processor for computing absolute blood 
volume using a reconstructed image acquired with the blood signal nulled by an inversion 

recovery preparation. 

FIGURE 5 shows a block diagram of a processor for computing a difference 
image having enhanced contrast for vascular space occupancy abnormalities. 

Detailed Description of the Preferred Emb odiments 

With reference to FIGURE 1, a magnetic resonance imaging scanner 10 
includes main magnet coils 12, which are preferably superconducting coils, although resistive 
main magnet coils or a permanent magnet can also be employed. The main magnet coils 12 are 
energized to generate a substantially uniform main magnetic field in an examination region 14. 
Magnetic field gradient coils 16 produce gradients in selected spatial directions to spatially 
encode magnetic resonances that are generated by energizing a radio frequency coil 18. In 
FIGURE 1, a whole-body radio frequency coil 18 is shown; however, local coils such as head 
coUs, phased radio frequency coil arrays, SENSE coils, and the like can be used instead of or in 
conjunction with the whole-body radio frequency coil 18 to excite magnetic resonances and/or 
to detect magnetic resonance echoes. 

A magnetic resonance sequence controller 30 coordinates and controls a radio 
frequency transmitter 34 that is coupled to the whole-body radio firequency coil 18 or another 
radio frequency coil to excite magnetic resonance echoes, and controls magnetic field gradient 
controllers 32 coupled to the gradient coils 16 to spatially encode the excited magnetic 
resonance echoes. One or more radio frequency receivers 36 coupled to the whole-body radio 
frequency coil 18 or another radio frequency coil detects, demodulates, and digitizes the 
magnetic resonance echoes and stores digital magnetic resonance samples in a k-space memory 



40. A reconstruction processor 44 peafonns a Fourier transform-based image reconstruction or 
other type of image reconstruction to generate one or more reconstructed images from the 
stored k-space magnetic resonance samples. 

The reconstructed images are stored in an image memory 46, processed by a 
video processor 50 and displayed on a user interface 52, transmitted over a local computer 
network or the Internet, or otherwise processed. Preferably, the user interface 52 includes a 
display, printer, or other output device that aUows a radiologist or other operator to view, 
render, or otherwise manipulate the reconstructed images. Moreover, the user interface 52 
preferably enables the radiologist or other operator to communicate with the magnetic 
resonance sequence controller 30 to create magnetic resonance imaging sequences, modify 
imaging sequences, execute imaging sequences, or otherwise control the magnetic resonance 

imaging scanner 10. 

With continuing reference to FIGURE 1, for clinical or diagnostic applications 
in which the blood volume is of interest, a measure of the vascular space occupancy (VASO) is 
perfonned. The approach uses an inversion recovery magnetic resonance excitation sequence 
having an inversion time 60 optimized for substantial nulling of the blood signal based on a 
blood Tl value particular to a certain magnetic field strength and hematocrit range. The 
inversion recovery magnetic resonance excitation sequence substantially nulls the signal from 
blood so that the magnetic resonance signal corresponds predominantly to tissue, with 
substantially negligible contribution of blood. This approach takes advantage of a substantial 
independence of the Tl value of blood on blood oxygenation, which allows the inversion 
recovery blood nulling to be operative for different types of vessels (arterial, arteriolar, 
capillary, venular, and venous) as well as throughout a physiological perturbation cycle. 

A suitable value for the mversion time (Tl) 60 for nulling blood can be 
obtamed m a number of ways. In one suitable approach, a blood Tl measuring sequence 62 is 
applied by the sequence controller 30 to measure a Tl value of a representative blood sample 
using a blood perfusion apparatus, and the inversion time (Tl) 60 is computed therefrom. 
Alternatively, a major blood vessel can be identified in a reconstructed image, and a direct 
measurement of the Tl value of blood inside the identified blood vessel is obtained using 
conventional magnetic resonance imaging sequences. In another approach, as it has been found 
in practice that the Tl of blood is generally not subject-dependent for most human subjects 
having substantially normal blood, suitable inversion times (Tl) can be computed from 
measurements of a representative blood sample of substantially normal human blood or of 
animal blood of suitable mammalian origin and collected in a table 64, such as Table I and 
Table U contained herein. These values are generally suitable for substantially normal human 
blood, and include inversion time dependence upon the main magnetic field strength (Table I 
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providing values for a 1.5T field, and Table H providing values for a 3.0T field) and on the 
sequence repeat time (TR). Rather than employing tabulated data, the inversion time 60 can be 
related to magnetic field strength, repeat time TR, and optionally other parameters by an 
empirical functional relationship or other suitable relationship. The tabulated TI values are 
preferably used as a guideline for determining an inversion time TI that substantially nulls 
blood so that signals are predominantly from tissue. The tabulated values do not exclude other 
TI values that can also accomplish blood nulling. 
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TR (ms) 


optimal TI (ms) 


remaining qray signal 


1000 


409.4499457 


0.039848653 


1500 


551.5765181 


0.07104825 


2000 


659.2454703 


0.100852361 


2500 


738.9402994 


0.126845436 


3000 


796,8489702 


0.148292976 


3500 


838.3345867 


0.165337186 


4000 


867.7421743 


0.178518564 


4500 


888.427932 


0.188505318 


5000 


902.8982691 


0.195951912 


5500 


912.9810254 


0.20143462 


6000 


919.987164 


0.205430439 


6500 


924.8460607 


0.208318508 


7000 


928.2112602 


0.21039169 


7500 


930.5397637 


0.21187148 


8000 


932.1498927 


^ 0.212922716 


8500 


933.262774 


0.213666531 


9000 


934.0317304 


0.214191053 


9500 


934.5629341 


0.214559874 



Table I - Optimal inversion times (TI) for blood at 1.5 Tesla magnetic field 
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I09000/00 


ouuu 
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0.195801064 


6500 


1101.689382 


0.20057404 


7000 


1109.605744 


0.204255213 


7500 


1115.458796 


0.207077192 


8000 


1119.781556 


0.209229423 


8500 


1122.971521 


0.210863669 


9000 


1125.324122 


0.212099926 


9500 


1127.058395 


0.21303208 



Table H - Optimal inversion times (TI) for blood at 3.0 Tesla magnetic field 



An inversion time adjustment is optionally performed to compensate for 
deviations in blood Tl value resulting from abnormal hematocrit values, sickle-cell 
pathologies, or another blood abnormality of a specific imaging subject. In one suitable 
adjustment technique, the blood-nulling inversion recovery imaging sequence is performed for 
several inversion times around the blood nulling inversion tune selected from Tables I and E, 
and the blood nulling inversion time 60 is selected as the inversion time providing substantially 
negligible image signal from a large blood vessel. A data flow line 66 in FIGURE 1 
corresponds to selecting the blood nulling inversion time 60 based on magnetic resonance 
measurements performed using the scanner 10. 

FIGURE 2 diagrammatically shows a suitable magnetic resonance imaging 
pulse sequence employing a blood-nulling inversion recovery magnetic resonance excitation 
sequence 70 to null the blood signal, and an exemplary single shot echo planar imaging readout 
72. Note that the echo planar imaging readout is exemplary only, and does not exclude 
employing additional or other magnetic resonance imaging, magnetic resonance spectroscopy 
or localized spectroscopy detection schemes. For example, the readout sequence 72 can be a 
single-shot imaging sequence, a single-shot echo planar sequence, a muUi-shot imaging 
sequence, a spectroscopy sequence, a multiple slice image, a one-dimensional, two- 
dimensional, or three dimensional spatial encoding sequence, a fractional k-space acquisition 
sequence, a spin echo readout sequence, a gradient echo readout sequence, or the like. 



Lx)oking first at the inversion recovery excitation sequence 70, an inversion 
pulse 74 is applied to invert the spins. Preferably, the inversion pulse 74 is a 180° pulse 
implementmg a 180° flip angle for the spins. In this case, the inversion time TI 60 is 
substantially representative to induce proper blood signal reduction based on the Tl value of 
blood for the particular experimental conditions (field, hematocrit, oxygenation, etc). It is also 
contemplated, however, to use an inversion pulse having a flip angle greater than 90° but other 
than 180°, in which case the appropriate inversion time TI 60 is readily computed from the 
proper Tl and acquisition parameters such as the repetition time TR. In general, the inversion 
tune TI 60 is selected as a time during which the longitudinal component of the flipped spins of 
blood decay from the flip angle to the crossover or null position. The null condition 
corresponds to a zero-crossing pomt of the longitudinal spm component; as the longitudinal 
spin component decays from the flipped or inverted alignment back toward the normal, 
non-inverted alignment it passes through a point where the longitudinal component passes 
substantially through zero, that is, the longitudinal spin component through a substantially zero 
crossing pomt. Substantially zero is understood to correspond to a substantially negligible 
blood signal such that the acquired magnetic resonance signal predominantly contams signals 
from tissues other than blood. 

Preferably, the inversion pulse 74 is not accompanied by a spatial encoding 
magnetic field gradient pulse or is accompanied by a relatively small spatial encoding magnetic 
field gradient pulse. This ensures that the spins of blood throughout the subject region of 
mterest reach the null condition after the inversion time delay TI 60. In particular, the blood 
nulling is independent of blood flow since the blood-nulling inversion pulse 74 is spatially 
non-selective or selects a relatively large region. Thus, flowing blood that flows into the slice 
of mterest at the time of excitation or at the time of readout is nulled appropriately. 

At the blood null condition, that is, after a delay time corresponding to the 
inversion time TI 60, an excitation pulse 80 is applied in conjunction with a slice-selective 
magnetic field gradient pulse 82 to excite spins in a selected slice of the subject region of 
interest. The excitation pulse 80 is preferably a 90° excitation pulse having a flip angle of 90°; 
however, an excitation pulse with other than a 90° flip angle is also contemplated. Because the 
blood is in a null condition at the time the excitation pulse 80 is applied, negligible magnetic 
resonance signal is excited in the nulled blood by the excitation pulse 80. Tissue such as fat, 
gray and white brain tissues, and the like generally have a different Tl value from that of 
blood, and so these tissues are not at a null condition at the time the excitation pulse 80 is 
apphed. Hence, the excitation pulse 80 excites magnetic resonance predominantly in tissue. 

The exemplary single shot echo planar imaging readout 72 samples the 
magnetic resonance excited in the tissue by the excitation pulse 80. The illustrated single shot 
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echo planar imaging readout 72 is a conventional readout including a 180° spin refocusing 
radio frequency pulse 84 and sUce-selective gradient pulse 86 that create a spin echo in a 
selected slice at a time-to-echo interval TE after the excitation pulse 80. A series of 
phase-encoding magnetic field gradient pulses 88 and a generally sinusoidal read magnetic 
field gradient waveform 90 step through a grid of k-space values in the selected slice while the 
radio frequency receiver 36 of FIGURE 1 performs samplmg 92 of the spin echo. The magnetic 
resonance signal k-space samples are stored m the k-space memory 40 of FIGURE 1 and are 
processed to produce a reconstructed unage representation as described previously. 

The single shot echo planar unaging readout 72 shown in FIGURE 2 is 
exemplary only. Substantially any type of magnetic resonance imaging or magnetic resonance 
spectroscopy readout or acquisition sequence can be employed. For imaging applications, the 
readout sequence may be a fast, single-shot sequence that acquires at least one slice per 
excitation. For example, the exemplary smgle shot echo planar imaging readout 72 can be 
replaced by a short gradient echo readout sequence. An echo time (TE) of the readout shown m 
FIGURE 2 is preferably kept short so that the readout is performed while the blood remains 
substantially in the nulled condition and to minimize contributions such as those from the 
blood-oxygen-level-dependent (BOLD) effect. Optionally, one or more additional inversion 
pulses are appUed to keep the blood close to the null condition during longer readout 
sequences. 

Although the inversion recovery blood-nulling sequence 70 is preferred, other 
pulse sequences can be employed to substantially null the blood signal so that the magnetic 
resonance is predominantly due to tissues other than blood. For example, the use of a long echo 
time (TE) at high magnetic field produces a substantially blood-nulled signal; however, in 
sequences of this type the extravascular blood oxygenation level dependence (BOLD) 
contribution is large, complicating determination of other physiological parameters such as 
cerebral blood volumes (CBV) and the like. In one suitable approach for overcoming this large 
BOLD effect, the CBV contribution is quantified from interpolation of the TE dependence 
from a series of TE values at sufficiently long TE to avoid intravascular contributions that 
would occur at shorter TE and/or from differences between such interpolations under different 
physiological conditions or between normal and diseased tissue. 

The blood-nulling magnetic resonance sequence of FIGURE 2 is optionally 
repeated with a repetition time (TR) to acquire reconstructed images for multiple sUces or for 
multiple three-dimensional (3D), two-dimensional (2D), or one-dimensional (ID) spatial 
encodings, for example to obtain a three-dimensional reconstructed image volume. In 
FIGURE 2, an exemplary succeeding spatially non-selective inversion recovery pulse 74' 
corresponds to the beginning of the next repetition of the blood-nulling imaging sequence. 



Optionally, because the blood nulling establishes a substantially negligible blood signal, it is 
possible to acquire multiple slices or multiple ID, 2D, or 3D spatial encodings per single 
nuUmg condition. This applies to both imaging and spectroscopy applications. 

Reconstructed images acquired using inversion-recovery blood-nulling 
magnetic resonance sequences such as the exemplary sequence of FIGURE 2 can be used in 
various ways. For example, while it is known in the art to process a reconstructed image to 
determine proton density (p), Tl, and T2 weightings, the determined weighting values are 
typically composite values having contributions from both blood and tissue. By applying 
substantially the same processing to the blood-nulled reconstructed images, more accurate p, 
Tl, and T2 weightings for tissue can be obtained without interference from the blood signal. 
Optionally, other tissues are nulled in addition to blood, to minimize contributions from those 
other tissues or components. For example, such combined blood and tissue nulling can be used 
to largely isolate a magnetic resonance signal from cerebral spinal fluid (CSF) in the brain. 

Reconstructed images acquired using inversion-recovery blood-nulling 
magnetic resonance sequences contain information pertaining to vascular space occupancy 
insofar as the images includes contributions from tissue but substantially exclude contributions 
from the blood volume. By applying a physiological perturbation, the blood volume can be 
changed. For example, cerebral blood volume undergoes vasodilation responsive to visual 
stimulation and breath-hold. Similarly, cerebral blood volume undergoes vasoconstriction 
responsive to hyperventilation. Blood volume changes can also be induced by administration of 
a selected drug, for example, but not limited to, for stress testing or the assessment of vascular 
compliance. Moreover, certain diseases such as cardiac ischemia, stroke, cancer, vascular 
deformations, and the like, represent a chronic or transient physiological perturbation that can 
cause a change m blood volume detectable with the methodology described herein. 

In most of the above cases, the change in blood volume occurs principally in 
the microvessels rather than in the large vessels and parenchyma. Advantageously, the blood 
volume changes measured using the blood-nulled reconstructed images reflect the parenchymal 
blood volume (denoted BV herein) which substantially coiresponds to the volume of the 
microvessels without contributions from the larger vessels. In contrast, vascular volume effects 
measured by BOLD and many other existing techniques include the large vessels and 
parenchyma and other tissues close to these vessels (for example, cerebral spinal fluid). For 
parenchymal tissue (tissue with perfused blood), the parenchymal vascular space occupancy 
(VASO, also denoted ^ herein) is given by: 
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=^ (1), 

BV + V V 

^ ^ ^ tissue par 

where BV is the blood volume, Vdssue is the pure tissue volume (without blood), and Vpa^ is the 
volume of the parenchymal tissue, that is, Vpar=BV+Vtissuc. The measured parenchymal vascular 
space occupancy ^ advantageously is more sensitive to physiological perturbation, including 
permanent disease-induced perturbation, than is the total vascular volume. Changes in 
large-vessel volume (that is, outside parenchymal regions) are also accessible by the blood 
nulling approach. At appropriate resolution, changes in large-vessel volume do not interfere 
with the indicated parenchymal blood volume changes due to the applied spatial encoding. 

For parenchymal tissue, the magnetic resonance signal (denoted by S herein) is 
proportional to a sum of the magnetization contributions of the microvessels and of the pure 
tissue. In general, the total signal Spar from the parenchymal tissue (including both tissue and 
perfused blood) can be written as: 

Spar = Si,,^^ + Sjissy^ (2)> 

where Sbiood and Sussue are the blood and tissue signals, respectively. Inserting the vascular space 
occupancy ^ and water density factors Cpar and Cbiood for parenchymal tissue and blood, 
respectively, into Equation (2) yields: 

where T2biood and T2tissue are suitable time constants for the exemplary spin echo readout 72 of 
FIGURE 2 due to the spin refocusing produced by the radio frequency pulse 84. For a gradient 
echo readout, the T2 values should be replaced by T2* values in Equation (3). The terms Mbiood 
and Mtissue correspond to the initial transverse magnetization of the blood and the tissue, 
respectively, produced by the excitation radio frequency pulse 80. As a consequence of using 
the blood-nulling excitation sequence 70, however, the term Mwood corresponding to the blood 
is substantially zero, and so Equation (3) reduces to: 

Sp„ s iC^ -4-C„„od)-M,3s« -e--^^- (4). 
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A fractional parenchymal signal change AS/S occurring responsive to a physiological 
perturbation, defined as: 



AS S^-S;^ 



'par 



(5). 



can be written by combining Equations (4) and (5) as: 



AS _ (C^ -C,^)-e-^^- ^(C,, -C,^)-e-^- ^^^^ 



where the superscript "act" references the perturbed state and the superscript "rest" references 
the unperturbed rest state or another suitable reference state. Dividing the numerator and 
denominator of Equation (6) by Cpe^-'^^'^''""''^''^ the fractional parenchymal signal change 
AS/S can be written as: 



(1 _ ^act ^ ^blood ^ . g-TE/(l/T2Sue-l/T2S«r) - (1 - . ^blood ) 

AS ^par ^Par (7)^ 

S ^2 ^rest ^ blood j 

Cpar 



Making use of the series expansion approximation e'^-l+x and defining an apparent change m 
transverse relaxation rate of tissue AR2 as: 



AR2 = 



' ' (8), 



tissue tissue 



Equation (7) can be rewritten as: 



-AR2.TE - 4^^^ .^blood ^^act . C blood .aR2>TE + 4^^ .:^blood 
AS ^par ^par ^par ^gy 



( 1 _ ^rest ^ ^ blood ^ 
^par 
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Defining a rate of change of the vascular space occupancy as A£,=^^''*-^"'' and neglecting 
term (Cwood/Cpar) AR2 TE in Equation (9) and simplifying yields: 



-AR2 TE - 

AS^ ^par (10). 

S ^ J _ ^rest ^ ^ blood ^ 

^par 



FIGURE 3 shows a block diagram of an exemplary processor 100 that 
computes a rest blood volume BV'' and a blood volume change ABV based on Equation (10). 
The processor 100 receives imagmg data input from the image memory 46. The imaging data 
includes images taken at a plurality of echo times TE 102 with the imaging subject at rest, and 
corresponding images taken at the plurality of (that is, two or more) echo times TE with the 
imaging subject in a perturbed state. A processor 104 computes the fractional signal difference 
AS/S between the unage of the subject in the rest state and the image of the subject in the 
perturbed state for each echo time TE to form a table 106 of fractional signal difference AS/S 
versus echo time TE. Recognizing that Equation (10) is a linear equation in echo time TE, a 
linear regression processor 110 computes a slope 112 and an ordinate-intercept 114 of the AS/S 
V. TE relationship. 

A rest blood volume BV'* calculator 116 computes the rest blood volume 
gyrest ^Yic slope 112 according to: 



slope = r 



-AR2 



l-# 



rest Cblood 



1- 



BV 



rest 



(11), 



^blood 



'par 



"par 



where Equation (1 1) is the slope component of the linear relationship between AS/S and TE .of 
Equation (10). The rightmost side of Equation (11) is written with the rest parenchymal 
vascular space occupancy ^"^'^ replaced by blood volume BV'' divided by the volume Vpar of 
the parenchymal tissue under consideration in accordance with Equation (1). The BV*^* 
calculator 116 uses literature values for the water density factors Cpar and Cbiood, the known 
volume Vpar of the parenchymal tissue under consideration, and a value for AR2 calculated 
from T2tissue,act and T2tissue,rest values obtained by fitting the signal value S at several echo times 
TE for each of the perturbed and rest states, respectively. If Vpar is unknown, the method can 
be used to determme the absolute blood volume fraction or vascular space occupancy ''^^ ). 
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A blood volume change ABV calculator 122 computes the blood volume rate 
of change ABV 124 from the ordinate-intercept 114 according to: 



V c 
Intercept = . ^ ^ 



v'~V C , 

V par J 

where Equation (12) is the ordinate-mtercept component of the linear relationship between 
AS/S and TE of Equation (10) written with the rest parenchymal vascular space occupancy 
replaced by BV"''A^par and the change in vascular space occupancy A^ is replaced by 
ABVA^par. The ABV calculator 122 uses literature values for the water density factors Cpar and 
Cbiood, the known volume Vpar of the parenchymal tissue under consideration, and the rest blood 
volume BV^ 118 computed by the BV"=' calculator 116. Alternatively, the absolute blood 
volume fraction or vascular space occupancy C^*^' ) can be used if Vpar is unknown. 

FIGURE 4 shows a block diagram of another processor 130 that computes an 
absolute blood volume BV. Computation of BV by the processor 130 is not based on a 
transient response to a physiological perturbation. Rather, data from the image memory 46, 
which may be a single image, is processed using Equation (4). The image is preferably 
acquired using a short echo time TE, so that the term exp(-TEyT2&sue) in Equation (4) is 
suitably approximated as unity so that the signal S is given by: 

par 

where the parenchymal vascular space occupancy ^ is replaced by the expression in Equation 
(1). The value of Mtissue at the inversion time TI 60 is related to longitudinal tissue relaxation 
time Tltissue according to: 



M tissue (TI) 



= 1 - 2e~'"^*«s^« + e'^'^'tissue (14), 



where TR is the repetition time of the pulse sequence (indicated, for example, in FIGURE 2). 
A suitable value for Tlussue can be measured using a known magnetic resonance technique such 
as saturation recovery pulse sequence and inversion recovery pulse sequence or, preferably, by 
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a series of inversion pulses so that blood remains substantially nulled, while tissue signal 
decays with the representative Tl^. In the case of brain imaging, it is known that the Tlussue 
values for white and gray brain matter differ. Hence, as shown in FIGURE 4 a TUhit. value 132 
and a Tlg„y value 134 are preferably determined. Notice that Sp„ also contains instrumental 
factors, for instance related to coil sensitivity and excitation homogeneity, the combined 
contribution of which we will call IF. This constant can also be determined during the Tl 
measurement from the voxels in ventricles, which only contains CSF and has a proton density 
very close to pure water. An M&sue processor 136 computes corresponding M^ssue values at the 
inversion time TI 60 for white and gray brain matter, and a normalization processor 140 
normalizes the image stored in the image memory 46 to produce a normalized image 142 
according to: 

S (tissue) = ^ (15)' 

where in Equation (15) the tissue magnetization M&sue is selected as the magnetization term 
corresponding to the Tlg„,y value 134 for gray matter, and is selected as the magnetization term 
corresponding to the TUhite value 136 for white matter. A blood volume processor 144 
estimates the absolute blood volume BV 146 or the absolute blood volume fraction or vascular 
space occupancy i^"^') by applying Equation (13) to the normalized image 142. For pure 
parenchyma, an additional approach to determine blood volume is to use the tissue Tl value to 
null the tissue and only look only at remaining blood. Using the same IF determined m the 
blood nulling experiment, the signal in such a tissue-nulled voxel at short TE is given by: 

from which the BV (Equation (1)) and/or vascular space occupancy can be determined for 
known Cwood- 

FIGURE 5 shows a processor 150 that computes a clinical image from 
blood-nulled magnetic resonance images. The processor 150 is suitable in diagnostic 
appUcations for detecting an abnormaUty such as vasodilation or vasoconstriction caused by 
certain diseases such as cardiac ischemia, stroke, cancer, and the like. The image memory 46 
stores an image of the imaging subject to be diagnosed. The image is compared with a 
reference image 152 by a difference processor 154 to generate a clinical image 156. 
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In one suitable approach, the difference processor 154 computes an absolute 
difference between the images, so that normal areas appear dark in the clinical image 156 
whereas abnormal areas appear brighter in the clinical image 156 (or vice versa) due to 
differences between the subject image and the reference image in the vicinity of the 
abnormality. Optionally, the difference processor 154 computes a signed difference with a 
constant intensity level offset. In this approach, regions of vasodilation and regions of 
vasocontraction have opposite intensity polarities respective to the constant intensity level. 

The reference image 152 can be obtained from various sources. For example, 
an image of a normal subject can be employed for the comparison. Alternatively, in the case of 
a bilaterally symmetric anatomical structure such as the brain, a contralateral image can be 
used. For example, a suspect right-side of the brain can be compared with a presxmied normal 
left-side of the brain, preferably after a suitable left-right transposition of the contralateral 
comparison image. Still further, in cases of anatomical structures that are substantially 
homogeneous in the absence of an abnormality, a suspect portion of the organ can be compared 
with presumed normal portion of the same organ. 

In the case of brain diagnoses, the normalized image 142 of FIGURE 4 is 
optionally employed in the processor 150. The normalized image 142 advantageously has 
suppressed contrast between white matter regions and gray matter regions of the brain image, 
which makes the vascular space occupancy contrast more visible. 

In some clinical applications, large vessels changes may occur in addition to 
microvascular changes. These can also be detected, and, at sufficient spatial resolution, can be 
separated from the microvascular changes. 

The invention has been described with reference to the preferred embodiments. 
Obviously, modifications and alterations will occur to others upon reading and understanding 
the preceding detailed description. It is intended that the invention be construed as including 
all such modifications and aherations msofar as they come withm the scope of the appended 
claims or the equivalents thereof. 



